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Properties of bottomonia (T, \b an d T') in the Quark-Gluon Plasma (QGP) are investigated 
by assessing inelastic reaction rates and their interplay with open-bottom states (6-quarks or B- 
mesons) and color-screening. The latter leads to vanishing quarkonium binding energies at suffi- 
ciently high temperatures (close to the dissolution point), which, in particular, renders standard 
gluo-dissociation, g + T — > 6 + 5, inefficient due to a substantial reduction in final-state phase space. 
This problem is overcome by invoking a "quasifree" destruction mechanism, g, q, q+T — > g, q, q+b+b, 
as previously introduced for charmonia. The pertinent reaction rates are implemented into a kinetic 
theory framework to evaluate the time evolution of bottomonia in heavy-ion reactions at RHIC 
and LHC within an expanding fireball model. While bottom quarks are assumed to be exclusively 
produced in primordial nucleon-nucleon collisions, their thermal relaxation times in the QGP, which 
importantly figure into T-formation rates, are estimated according to a recent Fokker-Planck treat- 
ment. Predictions for the centrality dependence of T production are given for upcoming experiments 
at RHIC and LHC. At both energies, T suppression turns out to be the prevalent effect. 



I. INTRODUCTION 



One of the primary objectives of ultrarelativistic 
heavy-ion collisions is to produce and study the Quark- 
Gluon Plasma (QGP), a deconfined state of strongly in- 
teracting matter which is believed to have existed a few 
microseconds after the Big Bang. Recent observations 
from the Rclativistic Heavy- Ion Collider (RHIC) indicate 
that, in a temperature regime T ~ 1-2 T c (T c ~ 180 MeV: 
(pseudo-) critical temperature Q), the QGP is not the 
originally expected weakly interacting gas of quarks and 
gluons. These observations could, in fact, be related 
to recent results from lattice Quantum Chromodynam- 
ics (1QCD) which suggest the existence of resonance- (or 
bound-state) like hadronic states up to temperatures of 
T ~ 2T C , both in the heavy (QQ) 0,1111 and light (qq) 
quark sector 0, For charmonium states this opens 
the exciting possibility that a substantial part of the fi- 
nal yield arises from (secondary) recombination of c and 
c quarks in the QGP H Q (or at the hadronization tran- 
sition 0, H 0, 0), as opposed to the originally sug- 
gested signature of J/ip suppression 



The effective- 
ness of recombination mechanisms in statistical and ki- 
netic models sensitively resides on a thermal equilibration 
of the charm- (c-) quark momentum distributions in the 
QGP 0E1- On the one hand, kinetic c-quark equilib- 
rium is not easily conceivable based on perturbative QCD 
(pQCD) cross sections with thermal light quarks and glu- 



ons 0,0,0,0- On the other hand, a recent study 
based on resonant rescattering via "D" -meson states in 
the QGP 01 found substantially smaller c-quark ther- 
malization times (by about a factor of ~ 3 compared 
to pQCD estimates), being comparable to or even below 
the expected QGP lifetime at RHIC. Thus, the question 
arises whether J/ip suppression, or rather the absence 
thereof, is an appropriate QGP signature at collider en- 
ergies. At this point, the prospect of measuring T states 
at RHIC (and LHC) becomes particularly valuable. Due 
to their large mass, 6-quarks are not expected to kinet- 
ically equilibrate. Together with their small production 
cross sections at RHIC, the situation for bottom appears 
to be reminiscent to the charm sector at SPS, in the sense 
that J/ip suppression has been verified as the dominant 
effect, with small contributions from regeneration |?lll0j|. 

The main objective of this article is to quantify predic- 
tions for bottomonium production at collider energies un- 
der simultaneous inclusion of both dissociation and for- 
mation reactions. Earlier analyses of bottomonium in 
heavy-ion reactions have essentially focused on suppres- 
sion mechanisms, using either the Debye-screening pic- 
ture (combined with formation time effects) |19l \2(i |2lJ , 
or gluo-dissociation processes jH , g + T — > b + b [24| . 
In the present article we will assess the role of back- 
ward reaction channels leading to bottomonium regener- 
ation (secondary production) via bb coalescence through- 
out the QGP phase, as dictated by the principle of de- 
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tailed balance. Importantly, we incorporate effects of 
color screening via a substantial reduction of quarko- 
nium binding energies. This, in turn, forces us to re- 
place the gluo-dissociation cross section, which becomes 
inefficient (and ill-defined) for small binding energies, 
£ b < Aqcd — T, by "quasifree" scattering of ther- 
mal partons off the b and b quarks in the bottomonium 
state, g, q, q + T — > g,q,q + b + b. Despite an appar- 
ent suppression by one power of the strong coupling con- 
stant, ag, in-medium reduced binding energies render 
quasifree dissociation dominant over gluo-absorption due 
to a much larger phase space [T(il ]. 

Our article is organized as follows. In Sec.|H]we give a 
brief overview of production cross sections for both open 
and hidden bottom states in hadronic collisions, which 
will form our baseline for heavy-ion reactions via a scal- 
ing by the number of primordial (binary) collisions. In 
Sec. lIIII we evaluate equilibrium properties of bottomonia 
in the QGP (complemented by a brief motivation from 
recent 1QCD results in App. EJ>; we compute inelastic 
cross sections (Sec. lIII A")l and pertinent dissociation rates 
(Sec. IIII Bfl of T's and excited states (\b and T') from 
scattering off light partons, illustrating limitations of the 
widely used gluo-dissociation and suitable improvements. 
In Sec. lIVI we employ a kinetic-theory framework to study 
the time evolution of bottomonia in heavy-ion collisions, 
characterized by an underlying rate equation introduced 
in Sec. IIV Al In addition to inelastic reaction rates, the 
latter requires further input consisting of (i) a space-time 
(temperature) evolution of a heavy-ion reaction modelled 
by a thermal fireball (App. [5J) and (ii) the equilibrium 
limit of T densities (App. [UJ. Sec. IIV Bl contains our 
main results for the time evolution of T in central Au- 
Au (Pb-Pb) collisions at RHIC (LHC), with a discussion 
of the uncertainties related to the T equilibrium abun- 
dances relegated to appendices [D] an d El The evolution 
and feeddown contributions from excited bottomonia are 
addressed in Sec. IIV 01 In Sec. [V] we present our predic- 
tions for the centrality dependence of T and T' produc- 
tion at RHIC and LHC. We summarize and conclude in 
Sec. IVII including an outlook on future improvements of 
our analysis. 



II. INITIAL PRODUCTION OF 
BOTTOM(ONIUM) 

The initial production of bb pairs and related fractions 
of bottomonium states in hard (primordial) nucleon- 
nucleon (N-N) collisions constitutes an essential input 
for our subsequent calculations. Throughout the paper 
we will assume the total number of bb pairs in the system 
to be determined by hard production alone, since, due to 
the large 6-quark mass, m& 3> T, secondary production 
is expected to be negligible [25| (even under LHC initial 
conditions as considered below, m^/T > 4 — 5). The 
number of bb pairs at a given impact parameter, b, of an 
A- A collision (A=Au and Pb for RHIC and LHC, respec- 



tively) can therefore be written as 



jy _ _ a pp^bb 
^-inelastic 



N coll (b) R 



V I 



(1) 



where 



];Mlc denotes the total inelastic proton-proton 
t 42 mb and 78 mb for RHIC and LHC, 



cross section 

respectively) |26(, and N co u(b) is the number of pri- 
mordial N-N collisions estimated within the Glauber 
model. The inclusive cross section, cx^p^bi for bb pro- 
duction in p-p collisions is taken from recent estimates 
of Rcfs. 24, 27j. We use the values o^^y, — 2.0 /ib at 
RHIC (^/sWn = 200 GeV) and a pp ^ bl = 160 ^b at LHC 
{\/snn — 5.5 TeV), which in the latter case includes a 
factor 0.8 to account for shadowing corrections. Finally, 
R_ v =0.52 (0.29) for RHIC (LHC) denotes the fraction of 
bb pairs within a given rapidity interval (Ay ~ 3.6 for 
two thermal fireballs) which we have estimated from dis- 
tributions calculated in perturbative QCD |28j . 

The primordial numbers of bottomonia are not 
very well known; following Ref. 0, [2^, we employ 
fpp^T =3.5 nb (152 nb) for RHIC (LHC, where we have 
again included a factor of 0.8 shadowing correction). 

The initial number to be used for the thermal evolu- 
tion has furthermore to account for inelastic processes in 
the pre-equilibrium stages. We approximate these by the 
so-called nuclear absorption, which, in principle, can be 
extracted from data in p-A collisions. Again, this infor- 
mation is rather scarce. We therefore adopt the standard 
approach employed for charmonia, which consists of a 
constant (high-energy) nucleon-bottomonium absorption 
cross section, crj^, coupled with a Glauber model for 
the nuclear overlap function [2^|. We note that, at col- 
lider energies, coherence effects may induce substantial 
deviations from the naive collision scaling [3(j , typically 
amounting to a significant reduction in absorption. We 
therefore use rather conservative values for the absorp- 
tion cross section of a~^. = 3.1 mb at RHIC (4.6 mb at 
LHC) [24j, which, for simplicity, is assumed to be identi- 
cal for all bottomonium states. 

Different values for cr^uc w ^ affect the results for the 
centrality dependence of bottomonium production in A- 
A collisions as presented in Sec. [V] On the one hand, 
the number of "directly" produced bottomonia, defined 
as the primordial number subjected to suppression only, 
Ndir = iVprim S nuc Sqgp , depends on the nuclear absorp- 
tion cross section through the nuclear suppression fac- 
tor, S nuc (schematically, S nuc w exp[— pjv c^uc L(b)] with 
an effective nuclear density, pn —0.14 fm -3 , and impact 
parameter dependent path length, L(b)). On the other 
hand, the number of regenerated bottomonia is practi- 
cally independent of 
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III. UPSILON PROPERTIES IN THE QGP 

In this section we will assess inelastic cross sections and 
dissociation rates of the T and its excited states in the 
QGP. Our main objective is a combined and consistent 
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FIG. 1: (Color online) Gluo-dissociation (with vacuum bind- 
ing energy: dashed line, or with in-medium reduced binding 
energy: dotted line) and quasifree (solid line) cross section 
overlap with a thermal gluon distribution at T — 300 MeV 
(dash-dotted line). 



treatment of color screening of the heavy-quark poten- 
tial and inelastic reaction channels. In analogy to 
the present situation for (ground state) charmonia, and 
motivated by 1QCD results, we assume the bottomonium 
masses to be constant. The hadronic phase will be ne- 
glected altogether throughout this paper. Whenever pos- 
sible, we will try to motivate necessary assumptions by 
findings of recent 1QCD calculations (see App. [X] for a 
brief survey containing more details). 



A. Dissociation Cross Sections 

The most commonly used process to evaluate the bot- 
tomonium (or charmonium) dissociation cross section, 
a Y lss (Y=T, Xb, Y')j on partons is the QCD analogue of 
photo-dissociation, g + Y — » 6 + 6 [22,|2i|. In Fig.^ the 
dashed line represents for the ground state T(ls) 
with vacuum binding energy, e™ c ~ 1.1 GeV (defined 
with respect to the BB threshold, 2ms), as a function of 
gluon energy w in the T-rest frame. It is peaked around 
u> ~ 1.5 Eb — 1-5 GeV, which for moderate plasma tem- 
peratures, T ~ 300 MeV, is significantly above average 
thermal gluon energies, uj = 3T; in other words, the 
maximum of the cross section has rather little overlap 
with the phase-space weighted thermal gluon distribu- 
tion, p 2 f 9 (oj;T) (dash-dotted line in Fig. [TJ. While an 
increase in T and a reduction in binding energy (as ex- 
pected from color-screening) initially improves the over- 
lap, a small value of Eb (toward/close to the dissolu- 
tion point) eventually leads to a very narrow energy in- 
terval where the cross section is active (dotted line in 
Fig. and thus renders gluo-dissociation increasingly 
inefficient. Formally this is due to the fact that for a 
loosely bound (large-size) 6-6 system, the absorption of a 
hard gluon essentially occurs on a single 6-quark, in which 
case energy-momentum conservation cannot be satisfied, 



and higher orders in the multipole expansion have to be 
included. This problem is even more severe for the less 
bound (excited) bottomonia. 

To implement effects of color-screening we follow the 
approach suggested in Ref . |l(J, |3l| in the context of char- 
monia: for small binding energies we approximate In- 
states as two comoving 6 and 6 quarks, which individually 
interact with surrounding thermal par tons via lowest- 
order (quasifree) pQCD scattering [33], g(q,q) + 6 — > 
s(<?7 q) + fr (and likewise for 6). The reaction is considered 
to lead to Y breakup, g(q, q) + T — * g(q, q) + b + b, if the 
residual binding energy can be overcome (and conserve 
overall four-momentum) . Note that in this approach also 
thermal (anti-) quarks contribute to bottomonium dis- 
sociation via i-channel gluon-exchange. The i-channel 
singularities are regularized by a gluon Debye-mass, 
fiD ~ gT, with as — 0.26 as fixed in our analysis for 
charmonia (note that the t-channel gluon-propagator 
compensates one power in as, which formally renders 
the quasifree cross section of the same order as gluo- 
dissociation). The result for the quasifree dissociation 
cross section of T(ls) with in-medium binding energy 
(solid line in Fig.^l shows a much improved overlap with 
thermal distributions due to its monotonous rise with 
gluon energy. The threshold for the result shown is due 
to a combination of binding energy and thermal parton 
masses; for massive partons, necessary to reproduce the 
QGP equation of state, the efficiency of gluo-dissociation 
is even more suppressed if Eb < frig ~ gT. In addi- 
tion, the quasifree cross section is theoretically better 
controlled for small values of the binding energy than 
gluo-dissociation since its main contribution arises from 
gluons with relatively large (thermal) energies, oj ~ 3T. 
Its applicability therefore also encompasses excited bot- 
tomonia. 



B. Lifetimes 

To evaluate bottomonium suppression in the plasma, 
the inelastic cross section, ay ss , is converted into a disso- 
ciation rate Ty (or inverse lifetime Ty 1 ) by convoluting 
it with thermal distributions, f g ,q, of quarks and gluons, 

f q ^ k ,T)v leX ap s {s) . (2) 

Here, s = (q + k) 2 denotes the total center-of-mass 
(CM) energy squared of the collision of a Y at rest 
(q = (my,0)) with a thermal on-shell parton of four- 
momentum k. We repeat that throughout this paper we 
use the vacuum masses for the bottomonia (for simplicity 
we do not distinguish the three spin states of the \b but 
use their average mass). 

Let us start by discussing the lifetimes arising from 
gluo-dissociation using constant (temperature indepen- 
dent) Y binding energies for two different values of the 
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FIG. 2: (Color online) Y lifetimes vs. temperature in a 
QGP based on the gluo-dissociation process for two (T- 
independent) values of the bottom-quark mass and perti- 
nent binding energies according to Eq. Q. Solid lines: T, 
dash-dotted lines: T', dashed lines: \b- The light- (dark- 
) colored set of curves corresponds to mj = 5.280 GeV 
(m b = 5.018 GeV). 
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FIG. 3: (Color online) The dependence of Y binding ener- 
gies (symbols) on gluon screening (Debye) mass according to 
solutions of the Schrodinger equation with a screened Cor- 
nell potential 13 J . The lines indicate fits used for the nu- 
merical implementation in our calculations. The identifica- 
tion with temperature is provided by the perturbative ex- 
pression for the Debye mass (see text), amounting to, e.g., 
^ D {T C ) ~ 325 MeV. 



open-bottom threshold (2mb), cf. Fig. For the larger 
mb = 5.280 GeV (light colored curves), correspond- 
ing to the lightest B-meson, the lifetimes are rather 
well-behaved at temperatures relevant for RHIC (T < 
400 MeV), i.e., rather large for the T (solid line, around 
10 fm/c and above, which is significantly longer than the 
expected QGP lifetime), much smaller for the \b (dashed 
line, 1 fm/c and below), and still slightly smaller for the 
less bound T' (dash-dotted line). However, for temper- 
atures in excess of ~500 MeV, the T' lifetime becomes 
larger than the one of the x&: a nrs t indication of arti- 
facts induced by the gluo-dissociation cross section for 
small binding energies as discussed in the previous sec- 
tion. This trend becomes more pronounced upon reduc- 
ing the bottom-quark mass to mi, — 5.018 GeV (dark- 
colored set of lines in Fig. implying e\ — 13 MeV 
and £§' = 145 MeV. In this scenario, the T' lifetime is 
roughly two orders of magnitude larger than the T and 
Xb lifetimes, which is obviously unrealistic and clearly in- 
dicates deficiencies of the gluo-dissociation cross section. 

As a more systematic way to include the effects of in- 
medium heavy-quark potentials in our calculations we 
make use of the results of Ref. [33| obtained from a 
Schrodinger equation using a Debye-screened Cornell po- 
tential. Assuming a Debye mass of the form suggested by 
pQCD, [io = gT (with g ~ 1.8 as used in the quasifree 
cross section), the T {\b, T') binding energy is found 
to be reduced to about -550 MeV (150-200 MeV) at 
T c , reaching zero at about 4.3 T c (-1.8 T c ), cf. Fig. 
These values are in reasonable agreement with solutions 
to a Schrodinger equation based on recent (unquenched) 
1QCD potentials (except close to T c ) |34(. Since the 
screening close to T c might well be less pronounced, our 
in-medium scenario should be considered as an upper 
limit of a reduction in binding energies. The pertinent 
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FIG. 4: (Color online) Y lifetimes vs. temperature in a 
QGP based on gluo-dissociation with an in-medium bottom- 
quark mass according to in-medium reduced binding energies 
as computed in Ref. [33| (dark-colored curves). Solid lines: 
T, dash-dotted lines: T', dashed lines: \b- The light colored 
curves are identical to Fig. 0] i.e., for constant binding energy 
with m b = 5.280 GeV. 



effects in connection with the gluo-dissociation cross sec- 
tion are illustrated by the dark-colored curves in Fig. 0] 
(the light-colored curves are identical to Fig. Whereas 
close to T c the ordering is still reasonable, the strong re- 
duction of Eb with increasing T rapidly induces a very 
unrealistic behavior in terms of a reduction in the dissoci- 
ation rate (increase in lifetime). As discussed above, this 
artifact is induced by a large decrease in available phase 
space for the gluo-dissociation cross section which ren- 
ders bottomonium dissociation a very inefficient reaction 
as the temperature increases. 

As indicated in the previous section the (presumably 



5 




0.6 

T [GeV] 

FIG. 5: (Color online) Y lifetimes computed from quasifree 
dissociation using constant (light set of curves) and in- 
medium (dark set of curves) binding energies for the T (solid 
line), T' (dash-dotted line) and the xt (dashed line). 



unphysical) decrease in the reaction rates with increas- 
ing T can be remedied by introducing the quasifree dis- 
sociation processes which are not limited by phase space 
and therefore apply to excited states as well. The re- 
sults are summarized in Fig. |5J where the correspond- 
ing Y lifetimes are plotted as a function of temperature. 
The light-colored lines represent vacuum binding ener- 
gies with an open-bottom mass of nib — 5.280 GeV. The 
resulting lifetimes are somewhat larger than the ones us- 
ing gluo-dissociation (light-colored lines in Fig. [5J , but 
exhibit the (from the size and binding energies) expected 
hierarchy, Ty < r Xb < tt, at all temperatures. The real 
virtue of the quasifree approach manifests itself in a (pre- 
sumably more realistic) calculation using in-medium (T- 
dependent) binding energies corresponding to the dark- 
colored lines. Here, smaller eb facilitate bottomonium 
breakup, and the respective dissociation rates are always 
larger than the ones using the vacuum values for eb- In 
comparison to the gluo-dissociation results, two further 
remarks are in order: (i) At the highest temperatures the 
Y lifetimes arc larger than the ones calculated with gluo- 
dissociation using vacuum binding energies. The main 
reason for this is that the strong coupling constant, as, 
inherent in the b-b Cornell potential necessary to repro- 
duce bottomonium spectroscopy (which also enters the 
dissociation cross sections), is substantially larger (~0.5) 
than in the quasifree calculation (~0.26). After all, it is 
the Debye-screening of the potential which is related to 
an effectively reduced as in the medium, (ii) At tempera- 
tures T < 300(200) MeV the gluo-dissociation rates for T 
(T', Xb) using in-medium binding energies are not much 
smaller than the ones obtained with quasifree destruc- 
tion. Therefore, both contributions to the rate should, in 
principle, be accounted for (at higher T, gluo-dissociation 
with in- medium eb is negligible). However, since we have 
fixed as in the charmonium sector using the quasifree 
process alone to reproduce J/ ip data in Pb-Pb collisions 
at the SPS, we will refrain from adding the two contri- 



butions. In addition, as we will see below, most of the 
bottomonium suppression at RHIC and LHC occurs at 
temperatures above ^250 MeV. 

Unless otherwise specified, in the following we will use 
the bottomonium lifetimes pertaining to the quasifree 
mechanism with in-medium binding energies as our de- 
fault scenario (we recall that this should provide an upper 
estimate of the screening effects). For comparison pur- 
poses with existing literature, we will also discuss results 
using the gluo-dissociation reaction alone with vacuum 
binding energies. 



IV. KINETIC THEORY APPROACH 

To calculate the time evolution and final abundances 
of bottomonia in a high-energy heavy-ion collision, we 
employ a kinetic rate-equation approach, as has recently 
been applied to charmonium production |6|, |7j . On the 
one hand, such an approach has the appealing feature 
that it makes explicit contact to equilibrium properties 
of quarkonia |jj , as discussed in the previous section. On 
the other hand, it does not provide an accurate treatment 
for processes occurring far from equilibrium, in which 
case explicit transport calculations are preferable |35ll36| . 



A. Rate Equation 

In the limit that bottom quarks (and/or hadronic 
open-bottom states which might survive up to T ~ 2T C or 
so) are in thermal equilibrium with the bulk of the system 
(light partons), the time dependence of the Y number, 
Ny, obeys a simplified rate equation given by |37l ] 



dNy 



-r r (n y - N Y q ) 



(3) 



where spatial (temperature) gradients have been ne- 
glected. Here, Ty denotes the V-dissociation rate dis- 
cussed in Sec. IIII 51 above and N Y q the equilibrium num- 
ber. The first term on the right-hand side of Eq. © 
describes the loss of bottomonia due to dissociation, 
whereas the second term takes into account regeneration 
according to the principle of detailed balance, 



Xi + Y «-» X 2 + b + b 



(4) 



where X\ and X 2 are light-quark or gluon states. Ob- 
viously, the backward channel amounts to a recombina- 
tion (or coalescence) of a b- and 6-quark within the QGP. 
From Eq. © it is clear that the regeneration contribution 
to the Y yield is negligible if either the equilibrium num- 
ber is small compared to the actual number, N Y q <C Ny 
at all times, or the dissociation/recombination processes 
shut off altogether, i.e., the reaction rate T-f is very small. 
Eq. is amenable to an explicit solution which reads 



Ny(t) = S QGP (t) 



dr 



r F (r)V° q (7 
Sqgp{i~) 



(5) 
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FIG. 6: (Color online) Time evolution of temperature within 
our expanding fireball model for a massive parton gas in cen- 
tral A-A collisions at RHIC (dashed line) and LHC (solid 
line). The use of massless partons (with Nf = 2.5) affects the 
evolution rather moderately. 



10 



L0 



10' 



. 1 1 1 1 1 1 
f i\ 






LHC : 
Central A-A 




■ — ■ U qiidlKS, — J.ZOU VJC V 


S. \ 
r v.> s 

IRHIC N ^ 
r Central A-A ( 


— b quarks in-medium 

— Known B hadrons ; 

— All B hadrons : 

— b quarks in-medium 

E 



0.2 



0.4 



0.6 

T [GeV] 



o.x 



FIG. 7: (Color online) Y(ls)-equilibrium abundances for cen- 
tral A-A collisions at RHIC (lower 3 curves) and LHC (upper 
2 curves), according to different scenarios for open-bottom 
states in the QGP (see text and App. imi . 



where 



SqgpW = exp 



- / drr T (r) 



(0) 



is the QGP suppression factor representing the probabil- 
ity that an Y, present at t = 0, has survived at time t, 
and Ny is the number of primordially produced Y's after 
pre-equilibrium (nuclear) absorption. 

The solution to the rate Eq. © requires the time de- 
pendence of the temperature, T{t). To this end, we de- 
scribe the space-time evolution of the system with a sim- 
ple thermal fireball expansion [3^]. It is constructed to 
reproduce the main features (timescales and transverse 
flow velocities) of hydrodynamical simulations of A-A 
collisions 39], and has been applied before to thermal 
dilepton |40j and charm(onium) |l8l l3l| production at 
SPS and RHIC energies, cf. App.lBlfor more details. The 
resulting temperature profiles are summarized in Fig. 
the pure QGP phase lasts for about 3 fm/c at RHIC 
with initial temperatures, To, close to 400 MeV, while its 
duration is -8 fm/c at LHC with T = T(r ) ~ 1 GeV. 

The final ingredient needed for the solution of the rate 
Eq. is the equilibrium number, N Y q , of Y states, 
governing their regeneration contribution throughout the 
evolution. As elaborated in App. [U] for a fixed number 
of bb pairs in the system, N^ 1 depends on the spectrum 
of open-bottom states in the QGP. We consider the fol- 
lowing possibilities: 

(i) 6-quarks with m& = 5.280 GeV, 

(ii) 6-quarks with in-medium mass rrib(T) following from 
the relation 



2m fe (T) =m Y +e%(T) 



(7) 



with constant (free) Y"-masses and binding energies ac- 
cording to Fig. |3 

(iii) B-hadrons as listed by the particle data group [2f3| . 



and 

(iv) open-bottom hadrons as extrapolated from known 
charmed hadrons. 

The general trend is that a higher degeneracy and a 
smaller mass of available open-bottom states favors b 
and b quarks to reside in the former and thus entails 
a smaller equilibrium abundance of bottomonia, as re- 
flected in Pig- H (cf. also App.EU- 

So far our discussion has been residing on the assump- 
tion that open-bottom states have momentum distribu- 
tions in thermal equilibrium with the bulk matter. How- 
ever, due to the rather short duration of a few fm/c of 
the QGP phase in URHIC's, and a 6-quark mass that ex- 
ceeds expected early temperatures by a large factor, this 
is presumably not a realistic assertion 0, • Following 
a suggestion made in Ref. [3l]], all our results displayed 
in the main text below therefore contain a schematic cor- 
rection to account for incomplete 6-quark equilibration, 
see App. IG 21 with the pertinent thermal relaxation times 
taken from Refs. 0, 0] which in turn provide a reason- 
able description of recent RHIC data for single-electron 
spectra (attributed to D- and B-meson decays) 0, ■ 
The associated numerical uncertainty is illustrated by 
considering the thermal-equilibrium limit in App. [EJ 



B. Time Evolution of T in Central A-A 

With all ingredients in place we proceed to evaluate the 
time evolution of the number of T states in central Au- 
Au (Pb-Pb) collisions by convoluting the rate equation 
© over the thermal fireball expansion. As a reference 
scenario commonly used in the literature, we first con- 
sider T interactions using the gluo-dissociation process 
with vacuum binding energies. We then turn to our de- 
fault scenario corresponding to the quasifree mechanism 
with in-medium binding energies. 
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FIG. 8: (Color online) T(ls) abundance as a function of 
time for central (6=1 fm) Au-Au collisions at RHIC using 
the gluo-dissociation process (without feeddown from excited 
states). The solid line includes T suppression and regener- 
ation and is almost indistinguishable from the dashed line 
corresponding to T suppression only. The dash-dotted line 
represents the T-equilibrium number which governs the re- 
generation term in Eq. JIJ. 



FIG. 9: (Color online) T(ls) abundance as a function of time 
for central (b = 1 fm) Pb-Pb collisions at LHC using the 
gluo-dissociation suppression mechanism (without feeddown 
from higher bottomonium states). The solid line, including 
both suppression and regeneration, indicates substantial sec- 
ondary T production when compared to the dashed line (no 
gain term) , facilitated by rather large T-equilibrium numbers 
(dash-dotted line) and reaction rates. 



1. Gluo-Dissociation 



As elaborated in Sec. IIII Bl in the gluo-dissociation 
scenario we are limited to using vacuum binding energies 
which we define relative to the hadronic open-bottom 
threshold. This amounts to equating the bottom-quark 
mass, mj, to the mass of the lightest bottom hadron, i.e., 
nif, = tub = 5.280 GeV. This rather large mass entails 
comparatively large T-equilibrium numbers (cf. dash- 
dotted lines in Fig. (RHIC) and ED (LHC)), which for 
both RHIC and LHC significantly exceed primordial pro- 
duction in the later stages of the evolution. 

Nevertheless, at RHIC the effect of regeneration is es- 
sentially irrelevant (full vs. dashed line in Fig. |8J) due 
to small reaction rates at the relevant temperatures of 
T < 400 MeV, cf. Fig.0 

The situation is different at LHC (Fig. EJJ where sec- 
ondary production is found to be at the 40% level 
(dashed vs. full line Fig. EjJ, due to a combination of 
a longer plasma lifetime (facilitating the build-up of a 
larger "equilibrium" limit due to thermal relaxation of b- 
quarks) as well as the higher initial temperatures where 
the inelastic reaction rates are much increased (cf. Fig.[2J). 
However, within the gluo-dissociation framework the im- 
portance of regeneration at LHC crucially hinges on the 
large T-equilibrium abundance induced by the large mass 
of the bottom quark, m& = 5.280 GeV, which in turn 
was used to consistently implement the vacuum binding 
energy of the T. This behavior will be modified if in- 
medium binding energies are employed, as we will see in 
the following subsection. 



2. Quasifree Dissociation 

As discussed above, we use the quasifree approach 
to inelastic reaction rates to implement in-medium re- 
duced T binding energies following the estimates within 
a screened heavy-quark (Cornell) potential of Ref. [33| . 
Imposing Eq. (JJJ further entails appreciably reduced 
bottom-quark masses with correspondingly smaller T- 
cquilibrium numbers (see full line in Fig. 0). This, in 
turn, leads to A^. q being well below the initial hard pro- 
duction at both RHIC and LHC. At the same time, the 
reaction rates are much increased over the free gluo- 
dissociation scenario of the previous section. 

At RHIC, the T-equilibrium number is so small that 
secondary production is virtually absent, cf. dashed vs. 
solid line in Fig. 1101 However, the increased inelas- 
tic reaction rates lead to an appreciable suppression of 
nearly 50% of the initial number (after nuclear absorp- 
tion). This rather dramatic effect is almost entirely due 
to the reduced T binding energy induced by a screening 
of the b-b potential in the QGP (the binding energy is 
already reduced to ~550 MeV at T c , rapidly dropping 
further before reaching zero at approximately 4.3 T c ). 
Thus, the magnitude of T suppression at RHIC appears 
to be a rather sensitive measure for Debye screening of 
the in-medium heavy-quark potential. The effects be- 
come stronger at LHC, cf. Fig.^] Despite a significant 
increase of the final T number due to regeneration in cen- 
tral collisions (solid vs. dashed line), the overall effect is 
an almost seven-fold suppression of the primordial yield 
(after nuclear absorption), which again is mostly facil- 
itated by an increased reaction rate originating from a 
reduced binding energy, especially in the early phases. 

We furthermore check the sensitivity of our LHC cal- 
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FIG. 10: (Color online) T(ls) abundances as a function of 
time for central (6=1 fm) Au-Au collisions at RHIC us- 
ing the quasifree suppression mechanism with in-medium re- 
duced binding energies (no feeddown from higher bottomonia 
included). The solid line (T suppression and regeneration) is 
similar to the dashed line (T suppression only). This is a con- 
sequence of the small T-equilibrium abundances (dash-dotted 
line) . 




time [fm/c] 



FIG. 11: (Color online) The same as Fig. ^| but for central 
(6=1 fm) Pb-Pb collisions at LHC. While the total yield 
(solid line) builds up a large fraction of regenerated T's (due 
to noticable equilibrium abundances (dash-dotted line)), the 
prevalent feature is a net suppression by about a factor of 7. 



dilation to the input of the charged-particlc multiplicity 
which is controlled by the total entropy in the fireball, 
cf. App. [gl Reducing S (and therefore dN^/dy) by a 
factor of 2 and keeping everything else the same amounts 
to smaller initial temperatures (by a factor ~ 2 1 / 3 ), 
shorter lifetimes of QGP (5.5 fm/c vs. 8.5 fm/c; the 
duration of the mixed phase is changed very little) and 
smaller volumes. The latter imply larger open-bottom 
densities at given T and thus an increase (by a factor 
of 2) of the open-bottom fugacity, 7t,, which in turn in- 
creases the equilibrium limit of the bottomonium num- 
ber by that same factor. On the other hand, the smaller 
temperatures and lifetimes of the evolution also inhibit 
the thermal relaxation of bottom quarks, which reduces 



FIG. 12: (Color online) Same as Fig. Illl but with a total en- 
tropy of the fireball reduced to half (5=20000), corresponding 
to a charged-particle multiplicity per unit rapidity of about 
~1600. 

the correction factor TZ of Eq. (|C7|) . The combina- 
tion of these two effects leads to an equilibrium num- 
ber which both in magnitude and time-dependence is 
quite similar to the evolution scenario with the larger 
dN ch /dy ~ 3200 (5=40500). The main net effect is thus 
the less efficient dissociation due to the smaller temper- 
atures in the evolution leading to a final suppression fac- 
tor (relative to the initial yield after nuclear absorption) 
of ~l/4 for dA^ c h/dy ~ 1600 as compared to ~l/7 for 
dAch/dy — 3200. Furthermore, the final yield is now 
dominated by the suppressed direct contribution. 

We also note that at both RHIC and LHC most of the 
suppression occurs in the first 2-3 fm/c, underlining that 
T production is an observable which is sensitive to early 
(QGP) phases of central heavy-ion reactions. 



3. Comparison to Charmonium 

It is worthwhile to confront the bottomonium results 
with the ones for charmonium of our previous analysis 
along similar lines [7lll0|: when going from SPS to RHIC 
energies, we have predicted a transition from a charmo- 
nium suppression-dominated regime at SPS to mostly re- 
generation at RHIC, with an approximately constant net 
suppression relative to initial production. The situation 
for bottomonium found here is quite different: in scenar- 
ios with color screening we expect the dominant mecha- 
nism to be T suppression both at RHIC and LHC, which 
increases with collision energy. This suggests the inter- 
esting possibility of observing an absence of suppression 
(or even enhancement over initial production |ll| ) for 
charmonium at LHC |llj | together with substantial sup- 
pression of bottomonium which would provide a rather 
unambiguous signature of charmonium regeneration in 
particular, and QGP formation (including screening) in 
general. A similar feature (albeit less pronounced) might 
also be present at RHIC if color-screening and associated 
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prompt T(ls) 


~ 51%, 


T(la) 


from Xb{lP) decays 


~ 27%, 


T(ls) 


from y&(2P) decays 


~ 10%, 


T(la) 


from T(2s) decays 


~ 11%, 


T(la) 


from T(3s) decays 


~ 1%. 



TABLE I: Decomposition of the inclusive T yield into prompt 
and feeddown contributions from excited bottomonia in pp 
collisions at ^/sjvjv = 39 GeV ^| • 

reduction in T binding energies become operative below 
temperatures of ~ 2T C : in central Au-Au a 50% suppres- 
sion of (the much more tightly bound) T (after nuclear 
absorption) should be compared to a 30% net suppression 
of J/-0's after inclusion of the (prevalent) regeneration. 

We emphasize again that the magnitude of the net 
T suppression is quite sensitive to the strength of color 
screening (affecting the bottomonium binding energies, 
and, in turn, the 6-quark mass according to m& = (my + 
£b)/2): at both RHIC and LHC the scenario without 
screening (Figs. |S] and [j]J leads to little net suppression 
(contrary to Figs. 1101 andll II where screening is included). 
At RHIC, this is due to the larger binding energies lead- 
ing to a factor of ~10 decrease in the reaction rates 
around T — 300 MeV, while at LHC it is due to larger b- 
quark masses implying larger Y equilibrium abundances 
which facilitate the regeneration in the early phases. On 
the other hand, uncertainties in the implementation of 
the screening mass (e.g., using \io = \/l + Nf j 6 gT with 
Nf—3 rather than the default \in = gT) are small (less 
than 15% for the final abundances). We recall that all 
considerations so far concern the exclusive T yields in 
central A- A collisions. 



C. Excited Bottomonia and Feeddown 

A complete description of inclusive T production re- 
quires to account for feeddown contributions from (late) 
decays of higher bottomonium states. Thus, one also 
needs to determine the (in-medium) production system- 
atics of these states. In addition, it is expected that some 
excited bottomonia may be measured directly, providing 
further model tests. According to Ref. [24| and previous 
experimental studies of bottomonium systems at Fermi- 
lab 0], the final T yield observed in pp collisions at 
y/SNN — 39 GeV is composed of up to 50% from feed- 
down, see Tab. [I] For simplicity, we will in the following 
neglect T(3s) states and will not distinguish between Xb 
states. We will also assume that the composition of the 
primordial T yield is energy- and system-size indepen- 
dent. 

The dissociation rates for T' and Xb are evaluated 
within the quasifree formalism including in-medium 
binding energies following the screening estimates of 
Ref. [33 (amounting to e B ~150-200 MeV at T c , reach- 
ing zero at about 1.8 T c ), cf. Figs. [3] and (recall from 
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FIG. 13: (Color online) QGP suppression of bottomonia 
in central (6=1 fm) Pb-Pb collisions at LHC, using the 
quasifree destruction process. Solid line: T, dashed line: T', 
dash-dotted line: xt>- 

Sec. IIII 51 that the gluo-dissociation approach cannot be 
applied to the small binding energies at hand and there- 
fore does not allow for a consistent treatment of feeddown 
effects). For illustration purposes, we display in Fig. IT31 
the pertinent suppression factors, Eq. (0, for central Pb- 
Pb collisions at LHC (neglecting the gain term in the rate 
equation ©). The T' and the \b exhibit similar suppres- 
sion (dashed and dash-dotted (xb) line, respectively), al- 
most a factor of 10 stronger than for T mesons (solid 
line), which is, of course, due to their smaller binding 
energies. 

The much stronger suppression of excited bottomonia 
(as well as their smaller equilibrium limits) compared to 
direct T's has important consequences for the inclusive 
T yields at both RHIC and LHC which will be shown 
in the following section. For central A- A collisions, the 
suppression in inclusive T production is about 30-40% 
more pronounced than for exclusive yields as displayed 
in Sec. El 



V. CENTRALITY DEPENDENCE 

The results presented in the previous sections were cal- 
culated for central collisions (b = 1 fm). In this section, 
we extend our results to different centrality classes by 
solving the rate Eq. for Au-Au (Pb-Pb) collisions for 
various impact parameters, b, (as was done before for 
charmonia in Ref. Q), and then plotting the yields at 
the end of each time evolution as a function of the num- 
ber of nucleon participants, N part (b). 

To solve the rate equations, we focus on the inelas- 
tic reaction rates obtained from the quasifree suppres- 
sion mechanism for &-quarks in the QGP with in-medium 
mass nib(T) (corresponding to the lower set of curves in 
Fig. [5J and include thermal off-equilibrium corrections 
for N^(t) (according to eq. (jC7|) 1. The inclusive T yields 
account for the feeddown from T' and Xb for which we 
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FIG. 14: (Color online) Centrality dependence of Nr/N co n 
at RHIC (Au-Au collisions at -^/sjvjv = 200 GeV) using the 
quasifree suppression mechanism including color-screening. 
Dashed line: suppressed primordial T. Dash-dotted line: re- 
generated T. Solid line: total T. 

solve the same rate Eq. (01 with pertinent reaction rates. 
The primordial numbers of bb pairs and bottomonia (the 
latter subjected to an identical but centrality-dependent 
nuclear absorption) are obtained by scaling our input 
numbers of Sec. [n] according to the number of binary 
collisions, N co n(b) (characteristic for hard processes). 

The centrality dependence of the fireball evolution 
(App. [BJ) is constructed assuming the total entropy at 
given impact parameter to be proportional to N par t(b), 
Stot(b) oc N pav t(b) (characteristic for soft processes, 
which for RHIC data is approximately satisfied within 
20% or so; for LHC the uncertainty is larger due to com- 
peting effects anticipated from gluon saturation and an 
increasing importance of a A^ co n-component 44]). The 
initial transverse radius (rrj) is estimated in cylindrical 
symmetry from the nuclear overlap function with con- 
stant formation time, tranverse acceleration and longitu- 
dinal velocity, cf. Eq. (|B1|) (for the purpose of extracting 
the temperature and volume evolution, effects of finite 
ellipticities in noncentral collisions are not significant). 

In the following, we first study the centrality depen- 
dence of the T yield at RHIC, then at LHC, and finally 
the T'/T ratio. 

A. RHIC 

Fig. 1 141 shows our results for T production in ^snn = 
200 GeV Au-Au as a function of participant number, nor- 
malized to the number of binary collisions. As expected 
from the dependencies in Fig. ^3 the T yield (solid line) 
exhibits a substantial increase in suppression going from 
peripheral to central collisions, with negligible regener- 
ation (dash-dotted line), i.e., the final T yield is exclu- 
sively due to primordial production subject to suppres- 
sion in the QGP (dashed line). This is primarily a con- 
sequence of the small number of bb pairs created in hard 



FIG. 15: (Color online) Same as Fig. El but for Pb-Pb colli- 
sions at LHC (ysjvjv = 5.5 TeV) with diV c h/di/| ccn trai — 3200. 



nucleon-nucleon scatterings, implying an T equilibrium 
number which is well below the actual number at any 
time (more so toward peripheral collisions). The preva- 
lent effect is thus a significant suppression of the inclusive 
T yield with centrality, originating from (i) suppression 
of direct T's due to (substantial) color screening, and (ii) 
suppression of (less bound) excited bottomonia with sub- 
sequent feeddown. Overall, our findings for T production 
at RHIC are rather reminiscent of the J/*f? suppression 
pattern at SPS as calculated in a similar approach (in 
line with NA50 data |if|). 



B. LHC 

In Fig. 1151 we present our predictions for the central- 
ity dependence of the T yield at LHC within the same 
approach as for RHIC. The longer lifetime of the QGP 
phase and, in particular, the higher temperatures reached 
at LHC render the suppression of primordial T (dashed 
line) even more pronounced. Regeneration (dash-dotted 
line), however, is no longer negligible, accounting for up 
to 2/3 of the final yield in central collisions. This should 
leave further traces in pt-dependencies, such as a soften- 
ing of the spectra and a sizable v-i. Nevertheless, the 
situation at LHC is qualitatively similar to the one at 
RHIC as far as the centrality dependence is concerned, 
i.e., the prevailing manifestation of QGP formation is T 
suppression at both RHIC and LHC. 

We also illustrate again the impact of a smaller 
charged-particle multiplicity in Pb-Pb at LHC, leaving 
all other quantities unchanged, see Fig. As found for 
the time dependence in central collisions (Sec. IIV B 2(1 . 
the smaller initial temperatures and QGP lifetimes lead 
to a factor of ^2 less suppression. Also note that the 
absolute number of regenerated T's hardly changes, im- 
plying that the final yield is now dominated by the pri- 
mordial component. 
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FIG. 16: (Color online) Same as Fig. ^1 but with half the 
total entropy in the fireball at each impact parameter, corre- 
sponding to dNch/dy ~ 1600 in central collisions. 

C. T'/T ratio 

We finally extract the T'/T ratio at LHC which was 
already implicit in the centrality dependence of the inclu- 
sive T yield shown in the previous section. The smaller 
T' binding energy obviously facilitates its suppression, 
and its larger mass additionally inhibits regeneration 
(due to smaller equilibrium abundances). Consequently 
we find a decreasing T'/T ratio, by about a factor of 2 
(3) when going from peripheral (p-p) to central collisions, 
cf. Fig. El (feeddown to the T' has been neglected; note 
that the limit for p-p amounts to 33%, or 17% upon in- 
clusion of the dilepton decay branching ratios). We note 
that even for central collisions, our calculations amount 
to a ratio of T'/T ~9% which lies well above its thermal- 
equilibrium value of -5% (at T = 180 MeV) that one 
would expect if bottomonium states were entirely created 
via statistical hadronization ^(|. The reason is that the 
number of T' reaches its terminal value at the beginning 
of the mixed phase where the pertinent reaction rate, 
Tx', becomes small, while the equilibrium limit drops by 
a factor of two toward the end of the mixed phase (as for 
the T in Fig. ITT|) . On the other hand (see also Fig. fTTfl . 
the terminal number of T essentially happens to coincide 
with the equilibrium value at the end of the mixed phase 
(even though it is frozen much earlier). Therefore the 
final T'/T ratio is roughly a factor two larger than the 
equilibrium value at the end of the mixed phase, which 
is the result displayed in Fig.E|f° r central collisions. 



VI. CONCLUSIONS 

At heavy-ion colliders (RHIC and LHC) hadrons con- 
taining bottom quarks are expected to become an impor- 
tant new probe of the produced (novel) forms of strongly 
interacting matter. In this article we have presented a 
kinetic theory framework to assess the evolution of bot- 



FIG. 17: (Color online) Centrality dependence of Ny / /N? 
(supplemented by the dilepton branching ratios of 
B W (T')~1.3% and B MM (T)~2.5%, respectively; feed- 
down on T' states is not included) at LHC for Pb-Pb 
collisions at ^snn = 5.5 TeV. 

tomonia in the QGP phases of high-energy heavy-ion re- 
actions, neglecting any reinteractions in hadronic matter 
(which is presumably a fair approximation for T, Xb and 
T' states). The underlying rate equation requires sev- 
eral inputs: (1) initial bottomonium abundances (which 
we estimated from collision-scaled p-p data), (2) space- 
time evolution of temperature (for which we employed 
an expanding fireball model), (3) inelastic reaction rates 
of bottomonia, and (4) bottomonium equilibrium abun- 
dances. 

For the inelastic reaction rates we have evaluated two 
mechanisms for dissociation due to thermal light quarks 
and gluons in the QGP: (i) gluo-absorption (the analogue 
of photo-dissociation) and (ii) "quasifree" scattering of 
both quarks and gluons off individual b and 6-quarks 
in the bound state. Whereas processes of type (i) are 
appropriate (and expected to be dominant) for deeply 
bound quarkonia, processes (ii) allow to treat the case 
of small binding energies (where gluo-dissociation leads 
to anomalously small reaction rates). Consequently, 
when introducing effects of color-screening, the associ- 
ated (temperature-dependent) reduced binding energies 
mandate the use of the quasifree processes, while for vac- 
uum binding energies gluo-dissociation is employed. 

The equilibrium limits of bottomonia, which govern 
the gain term (regeneration) in the rate equation, follow 
the trends previously established for charmonia: for a 
fixed number of bb pairs in the system (at given central- 
ity), larger masses (and smaller degeneracies) of open- 
bottom states imply larger T-equilibrium limits (render- 
ing bottomonium formation more favorable). However, 
in most cases considered (in particular when implement- 
ing medium effects with relatively small 6-quark masses) 
the equilibrium limit was significantly below the initial 
value from hard production (after nuclear absorption), 
especially after applying a schematic correction for in- 
complete 6-quark thermalization. 
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Our main results may therefore be summarized as fol- 
lows: (I) color-screening (i.e., in- medium reduced bind- 
ing energies) leads to a substantial increase of the inelas- 
tic bottomonium reaction rates; (II) at RHIC, this man- 
ifests itself as a 50% suppression effect (factor ~2) in the 
QGP phase (with virtually no regeneration), whereas for 
vacuum T states almost no (thermal) suppression occurs 
(in which case suppression in inclusive yields is due to 
feeddown from excited states); (III) at LHC, assuming 
dN ch /dy ~ 3200 in central Pb-Pb, the net QGP sup- 
pression becomes substantially stronger (factor ~7 (1.4) 
for in- medium (vacuum) binding energies), with the final 
yield containing a 75% (30%) component of regenerated 
T's. If the charged-particlc multiplicity turns out to be 
smaller (e.g., dN c h/dy ~ 1600), the reduced QGP initial 
temperature and lifetime will lead to weaker net suppres- 
sion (e.g., factor ~2 less). 

Our general finding is thus that in heavy-ion reactions 
at both RHIC and LHC, T suppression is the prevalent 
effect as a function of collision centrality, and there is 
large sensitivity to color-screening in the QGP. This is 
quite different from our previous studies of charmonia 
where J/vP regeneration is expected to be the dominant 
mechanism for RHIC energies and above. 

Concerning future improvements, it is desirable to pro- 
vide a more accurate treatment of in-medium bottomo- 
nium properties, implementing both screening and disso- 
ciation processes, as well as non-thermal fo-quark distri- 
butions, simultaneously in a more microscopic approach. 
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APPENDIX A: BRIEF REVIEW OF LATTICE 
QCD RESULTS 

Finite-temperature 1QCD calculations exhibit a sub- 
stantial reduction of the heavy-quark free energy [^lisj . 

F Q Q(r,T) = U Q Q(r,T)-TS Q Q(r,T) , (Al) 

with increasing temperature. When implementing the 
free energy into a Schrodinger equation, J/ip (-0') mesons 
have been found to dissolve slightly above (below) T c • 
More recently, it has been argued that the internal 



energy, Uqq, is the more appropriate quantity to be 
identified with a potential. Corresponding solutions 
to the Schrodinger equation entail significantly larger 
dissociation tem pera tures of the J/ip of around Tdi SS — 
2T C HI H3, EJ III- This led, in fact, to better con- 
sistency with (quenched) lQCD results for charmonium 
spectral functions, in which the peak structures associ- 
ated with the lowest-lying resonances (J/ip and r) c ) dis- 
appear slightly below 2T C 0, 0, • With the mass (peak 
position) of the spectral functions approximately stay- 
ing constant |53|| . the key property of the charmonium 
spectral functions is thus a change in their widths, i.e., 
inelastic reaction rates. This is further corroborated by 
1QCD studies which find the J/ip width to increase across 
the phase transition, from very small below to about 
200 MeVjust above T c [5j| . Such a behavior is quite con- 
sistent |37jj with both hadronic model calculations below 
T c , which typically result in small inelastic reaction rates 
of <0.1 fm -1 (mostly due to ir,p + J/ip breakup) [54j . 
and parton- induced destruction above T c , typically of the 
order of 1 fin -1 . Also note that in the presence of color 
screening, which induces a reduction in the c-c binding 
energy (eg > 0), a constant charmonium mass implies a 
decrease in the open-charm threshold, according to the 
relation 

m# = 2m c — £b ■ (A2) 

A decrease in the open-charm threshold with T has in- 
deed been found in the 1QCD free energies (or poten- 
tials), at least above T c |55j . and complies naturally with 
an increase in the bound state width, due to larger phase 
space in the break-up reactions, e.g., g + J/tp — > c + c. 

Effects of the hadronic phase are neglected in this ar- 
ticle; first lQCD results (56J (as well as potential models 
based on lQCD |50j) indeed indicate much higher dis- 
sociation temperatures of the T than for J/ip, around 



APPENDIX B: THERMAL FIREBALL 
PARAMETERIZATION 

The starting point for the thermal fireball model is a 
cylindrical ansatz for the (proper) time dependence of 
the (isotropic) three-volume according to 

V FB (t) = (z + v z t)ir(r + ^-t 2 y . (Bl) 

The longitudinal expansion is characterized by (i) an ex- 
pansion velocity of v z — 1.4c representing a rapidity 
range of Ay ~ 1.8 typical for a thermal fireball width, 
and (ii) an initial length z = 0.6 (0.125) fm for RHIC 
(LHC) which carries the meaning of a formation (ther- 
malization) time, tq — Zq/ Ay =0.33 (0.07) fm/c at RHIC 
(LHC). The initial transverse size is determined by the 
centrality of the collision, e.g., r = 6.4 (6.6) fm for 
b = 1 fm Au-Au (Pb-Pb) collisions at RHIC (LHC). 
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If one further assumes the total entropy, S, of the sys- 
tem to be conserved (i.e., isentropic expansion), the tem- 
perature evolution of the fireball matter can be inferred 
by equating the entropy density to the one of a thermal 
medium. For a non-interacting QGP one has 



s(t) 



S 



V FB (t) 



d 3 k 



{l±f i (w i )]n[fi("i)} 



(2tt) 3 

[l T /*M]ln[lT/*M]} , 



(B2) 



where /j denote the Fermi-Dirac (i=q, q, upper signs) 
or Bose-Einstein (i=g, lower signs) distribution func- 
tions of (massive) quarks or gluons with on-shell ener- 
gies uii — \Jk 2 + mf (m,: thermal mass), and di stand 
for the pertinent spin-color-flavor degeneracies. We will 
consider two scenarios for the QGP equation of state: (i) 
massless partons with an effective number of quark fla- 
vors N f = 2.5; in this case, s = (d g + 10.57V/)47r 2 /90T 3 , 
and the temperature can be obtained explicitly as T(t) = 
const x (S/V FB (t)) 1/3 ; (ii) a three-flavor QGP with ther- 
mal parton masses, 



g 2 T 2 
6 



g 2 T 2 
6 



-g 2 T 2 , 



(B3) 

where mo =150 MeV is taken as a bare strange-quark 
mass; in this scenario, Eq. (|B2|I provides an implicit re- 
lation for T(t) which is solved numerically. However, we 
observe no significant differences for the temperature pro- 
files when calculating the entropy density for a QGP with 
massless degrees of freedom (and Nf = 2.5) or massive 
partons according to Eq. 1|B3(I . 

If the entropy density, s(t), drops below the critical 
value of the QGP phase, s°- GP = s QGP (T c = 180 MeV), 
we perform the standard mixed-phase construction, 



S 



= f 3 ^ + (l-f) s 



v FB {t) 



QGP 



(B4) 



with volume fractions / and 1—/ of matter in the hadron 
gas (HG) and QGP, respectively. The entropy density of 
the HG at T c , s HG (T c ) ~ 6.1 fm~ 3 , is computed from a 
HG including all resonances up to masses of about 2 GeV, 
with a baryonic chemical potential /j,b — 30 (2) MeV 
corresponding to (expected) chemical freezeout at RHIC 
(LHC). The charged-particle multiplicity is then deter- 
mined by fixing the total entropy of the system |38j. e.g., 
5=10000 (40500) for central collisions at RHIC (LHC) 
translating into dN^/dy ~ 800 (3200; we have evaluated 
a LHC scenario with 5=20000, i.e., dN ch /dy ~ 1600, 
which is closer to extrapolations based on present RHIC 
data [ill H3). As mentioned above, we neglect the 
hadronic evolution of the system. 



APPENDIX C: EQUILIBRIUM ABUNDANCES 

Following our assumption of a fixed number of b-b pairs 
(at given centrality) in the fireball, the V-equilibrium 
number not only depends on temperature and volume, 
but also on a bottom-quark fugacity tj, = 7^ via 



N?{T) = d Ylb (T) 2 V FB {t) 



d 3 p ( E Y (p) 

3 eX P 



(27T) 



T 



r (CI) 

(Ey(p) = y m Y + p 2 , dy- spin degeneracy). The tem- 
perature dependence of the fugacity is determined by re- 
quiring that at each moment the total number of b- (b-) 
quarks in open and hidden bottom states matches the 
primordially produced abundance 

l io(7b^ V op) 

where Ij are modified Bessel functions. The numbers of 
hadronic open and hidden bottom states are given by 

^op = V FB Y,n x , xe{B,B,B*,B*,A b ,...}(C3) 



N hid = V FB )^n y , ye{T,Xbo,Xbi,---}, (C4) 



v 



respectively, with 



a* P e ^(_^M) (C5) 



(27T) 



T 



the pertinent densities in Boltzmann approximation 
(d x . y : spin-isospin degeneracy). If the open-bottom 
states in the QGP consist of b- and &-quarks, Eq. 1C3(I is 
replaced by 



N 



op 



x e {b,b} 



(C6) 



As emphasized in Ref. [ij, for a small number of total 
bb pairs (N bb - < 1), it is quantitatively important in the 
determination of 7& to enforce exact bottom number con- 
servation (i.e., in each configuration) using the canonical 
ensemble, which is implemented in Eq. I|C2() by the ratio 
of Bessel functions. For large values of N bb , the ratio of 
Bessel functions approaches 1 , and one recovers the grand 
canonical limit. In Eq. IIC4|) . we include all bottomonium 
states listed in Ref. [2(|, but due to their relatively large 
masses compared to B-mesons (or bottom quarks), they 
play no significant role in the determination of 7^ via 
Eq. ljC2fl . The main uncertainty in this procedure resides 



in the number, N op , of open-bottom states, Eq. (|C3|) or 
(IC6|l , which depends on both their masses and degenera- 
cies, as elaborated in the following subsection. Another, 
more indirect uncertainty in N Y q (T) pertains to the as- 
sumption of thermal equilibrium of open-bottom states 
which is underlying Eqs. (|C1(1 through l|C6|l . A perti- 
nent (schematic) correction will be discussed in App. IC 21 
below. 
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1. Open-Bottom States in the QGP 

To assess the uncertainties associated with the open- 
bottom spectrum in the QGP we will distinguish two 
basic scenarios (with one further subdivision each). The 
first one is that of individual 6-quarks (with variations in 
their mass), as to be expected for a weakly interacting 
QGP. The second one is closer in spirit to a strongly inter- 
acting QGP, where (colorless) hadronic correlations per- 
sist above T c . Such a scenario has received renewed at- 
tention stimulated by findings of hadronic bound and/or 
resonance states in 1QCD simulations 0, Q (and effective 
models [HIHIIl) for T ~ 1-2 T c . Also, the strong col- 
lective phenomena (radial and elliptic flow) found exper- 
imentally at RHIC are not easily reconciled within a per- 
turbative framework. On the one hand, we will therefore 
investigate the survival of hadronic open-bottom states 
(most importantly B-mesons) with a varying degree of 
level density toward higher masses, i.e., either B-meson 
(and A;, baryons) as listed by the particle data group j2|| 
or additional states estimated from analogy to open- 
charm hadrons (since the expected initial temperatures 
at LHC are well above 2T C , we consider the "hadronic" 
scenarios only for RHIC). On the other hand, even if 
no hadronic states are formed in the QGP, at moder- 
ate temperatures T < 2T C , the strong coupling constant 
a s may not be small and induce significant temperature- 
dependences in the in-medium bottom-quark mass due to 
correlations with surrounding light partons. The main 
point in the following is that, at fixed N bb , differences 
in the open-bottom spectrum affect the 6-quark fugacity 
and thus, in turn, modify the V-equilibrium abundances. 
The general mechanism can be summarized as follows [jj : 
a larger number of available open-bottom states (either 
due to a large degeneracy or a small mass which fa- 
cilitates their thermal population) implies a smaller 7t, 
and thus a smaller equilibrium number of bottomonium 
states, A y q . In other words, heavier (and less) open- 
bottom states favor 6-6 quarks to reside in bottomonia. 

For definitcness, we consider the following four cases 
for open-bottom states in the QGP, summarized in Fig.0 

(i) b and b quarks with a fixed mass of m;, = 
5.280 GeV, corresponding to the mass of the light- 
est B-meson (LHC only, dash-dotted line in Fig.0; 

(ii) b and b quarks with in- medium masses mt,(T), 
where the temperature dependence is inferred in 
analogy to Eq. l|A"2"jl as 2m b (T) = m Y + £g(T) 
(Eq. 0) with the T-dependence of the T-binding 
energy, e£(T), taken from Ref. H^. As men- 
tioned above, here and throughout our work we 
assume all Y masses to be constant. The result- 
ing T-equilibrium abundances are indicated by the 
solid (LHC) and dash-double-dotted line (RHIC) 
in Fig. □ 

(iii) open-bottom hadrons as listed in the particle data 
book [2g: Ag, B, B* , B° s , Bf (RHIC only, dotted 



line in Fig. 01 . 

(iv) open-bottom hadrons including states with esti- 
mated masses as extrapolated from systematics of 
analogous (known) charmed hadrons (RHIC only, 
dashed line in Fig. . 

For the two 6-quark scenarios (i) and (ii) under LHC 
conditions (upper two curves in Fig.0, the main differ- 
ence amounts to a larger T equilibrium abundance for 
the larger 6-quark mass (case (i)); at high temperatures, 
where for case (ii) the 6-quark mass is essentially half 
the T mass (zero binding energy), the difference is less 
pronounced due to the flatter thermal distribution func- 
tions. Toward smaller temperatures (especially close to 
T c ), 7b increases significantly due to a decreasing JV* , 
and one notices the onset of the effect that, in the canon- 
ical ensemble (where b and b have to appear together), 
it becomes increasingly favorable to have a bb pair in an 
T state which has a lower mass than separate b and 6 
quarks (the vertical decrease of at T c is a conse- 
quence of the volume expansion during the mixed phase 
at constant temperature, implying a decrease in 7;,). 

At RHIC the 6-quark scenarios show a behavior for N^ 1 
analogous to LHC but with reduced magnitude mainly 
due to the reduction in N bb -. Close to T c , the scenar- 
ios with hadronic open-bottom states, (iii) and (iv), lead 
to a substantial increase in over the one with in- 
medium 6-quark mass, since the latter is well below 
the lowest B-meson mass. The sensitivity to higher- 
mass bottom hadrons becomes more pronounced with 
increasing temperature where they are more easily ex- 
cited, thus reducing the number of bottomonium states. 
This, in turn, implies that the 6-quark scenario (i), with 
m& = 5.280 GeV = raj (not shown in Fig0, is rather 
similar to the hadronic scenarios close to T c . 

We finally note that, during the mixed phase (i.e., at 
T e ), the open-bottom spectrum within each scenario is as- 
sumed not to change. With constant temperature and in- 
creasing 3- volume, Vfb(£), this implies that the bottom- 
quark fugacity, 7&(T, Vfb), decreases during the mixed 
phase (since N bb is conserved), and with it JVy q (T, Vfb)- 
The discontinuity in the first derivative of the tempera- 
ture evolution (cf. Fig. [fjj at the beginning of the mixed 
phase is thus reflected in %(T, Vfb) and N^(T, Vfb)- 

2. Incomplete 6-Quark Thermalization 

Due to the rather short duration of a few fm/c of 
the high-temperature (QGP) phases in URHIC's, and 
a 6-quark mass that exceeds expected early tempera- 
tures by a large factor, 6-quark thcrmalizatoin may not 
be a realistic assertion. Even for charm quarks (an ap- 
proach to) thermalization appears to be unlikely under 
RHIC conditions if only perturbative rescattering pro- 
cesses in the QGP are considered [l^. However, recent 
RHIC data j4j on the elliptic flow of (semileptonic) de- 
cay electrons which are attributed to D-mesons suggest 
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that charm quarks undergo substantial rescattering |13( 
in semi-central Au-Au collisions at y/s — 200 AGeV. A 
possible microscopic mechanism has been identified in 
Refs. |l7tll8j in terms of resonant c-quark scattering with 
light antiquarks via D-meson states surviving in the QGP 
at moderate temperatures, T< 2T C . It was found that, 
compared to perturbative calculations, the thermal re- 
laxation times could be reduced by up to a factor of 3-4, 
thus becoming comparable to the duration of the QGP 
phase at RHIC. A similar reduction was found for bottom 
quarks, but with the larger mass implying larger absolute 
values for the relaxation time scale by a factor of ~2-3, so 
that &-quark thcrmalization remains unlikely at RHIC. 

Consequently, the Y-equilibrium abundances shown 
in Fig. [7| should be improved by accounting for incom- 
plete thermalization. As in Ref. [3l|], we here follow 
a schematic relaxation time approach, by implementing 
into the thermal equilibrium abundances a reduction fac- 
tor 

K = l-exp(-f , (C7) 

where r cq is the bottom-quark equilibration time as ob- 
tained within the resonance-scattering model 17]. Typ- 
ical values of r oq range from r cq ~ 12 fm/c at T — 
300 MeV to r eq — 1.1 fm/c at T = 1 GeV, leading to a 
substantial decrease of the regeneration term in Eq. J2Jl, 
especially at RHIC. We will illustrate the uncertainties 
induced by this correction by considering the thermal 
limit in Add. IE! below. 



APPENDIX D: UNCERTAINTIES DUE TO 
OPEN-BOTTOM STATES 

The dependence of the Y-equilibrium numbers (with- 
out thermal relaxation correction) on different assump- 
tions about the open-bottom spectrum in the QGP was 
illustrated in Fig. [7| for the following scenarios: (i) b- 
quarks with m& = 5.280 GeV, (ii) &-quarks with in- 
medium mass rrih(T) , (iii) B-hadrons as listed by the par- 
ticle data group 26] , and (iv) open-bottom hadrons as ex- 
trapolated from known charmed hadrons. In this section 
we study the influence of scenarios (iii) and (iv) on the re- 
generation of T(ls), restricting ourselves to RHIC energy 
as it seems unlikely that hadronic correlations would per- 
sist at the temperatures anticipated at LHC (T < 1 GeV) . 
We will employ the rates corresponding to quasifree de- 
struction with in-medium binding energies (even though 
this implies an inconsistency between the rates and the T 
equilibrium number, it will allow us to estimate the max- 
imal possible effect as the quasifree rates are the largest 
considered in this work) . Fig. ^| pertains to the case 
where the open-bottom states in the QGP are given by 
the known £?-hadron spectrum. Since their mass tends 
to be larger than the in-medium 6-quark masses (as used 
in the main text, cf. Fig. I10|) . this scenario leads to an 
increase in the T-equilibrium abundances (light vs. dark 



- 2 

>• 

2 



o, 



RHIC 

7 Central Au-Au 
- Quasifree 


— Total '. 

— Direct 

. _ . IN 

eq 

Total - Known B hadrons 
. _ N - Known B hadrons - 

eq 
✓~ . 

/' " ' -- • ^ _ 


/ 

/' 

/' 


/ 


- , s ' 





) i ^ . -^ ., - i .-,.r,- \ --r-r; □ , , r-r 

1 2 3 4 5 6 



time [fm/c] 

FIG. 18: (Color online) Time evolution of T(ls) in cen- 
tral (6=1 fm) Au-Au collisions at RHIC assuming differ- 
ent open-bottom spectra in the QGP. Light curves: open- 
bottom states constituting known bottom hadrons with the 
dash-dotted curve showing the T-equilibrium abundance and 
the solid curve indicating the total T yield including regenera- 
tion. Dark curves: open-bottom states = in-medium 6-quarks 
as in Fig. ED 
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FIG. 19: (Color online) Same as Fig.rrglbut assuming a richer 
(hadronic) open-bottom spectrum in the QGP as estimated 
from known open-charm hadrons in the vacuum. 



dash-dotted curves in Fig. I18|) . However, the overall ef- 
fect of a slightly enhanced regeneration is small (upper 
solid line). 

If we assume the open-bottom spectrum to follow the 
same pattern as all known charmed hadrons, the increase 
in the T-equilibrium number over the in-medium 6-quark 
scenario is less pronounced (due to a smaller fugacity, 
7b, relative to using only known B-states), and, conse- 
quently, also the regeneration, leaving essentially no no- 
ticeable effect, cf. Fig. 1191 

Thus, overall, the effects of different open-bottom spec- 
tra in the QGP on T production at RHIC are very small. 
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FIG. 20: (Color online) Thermal off-equilibrium effects at 
RHIC for the suppression of T(ls) in central (6 = 1 fm) 
Au-Au collisions, using the quasifree process. The light set 
of curves corresponds to a full thermalization scenario with 
increased equilibrium abundances (upper dash-dotted line), 
leading to a 15% component from secondary production (up- 
per solid line). The darker lines correspond to Fig. HOl in the 
main text. 



APPENDIX E: UNCERTAINTIES DUE TO 
INCOMPLETE 6-QUARK THERMALIZATION 

Since 6-quarks (or -mesons) are not expected to kinet- 
ically equilibrate in central collisions at RHIC (and pos- 
sibly neither at LHC) 0, all previous figures included 
a thermal relaxation-time correction, Eq. (|C7ll . to the 
equilibrium limit of T states in the system (Figs. I5TH?|| . 
Due to the schematic character of its implementation, 
we study in this section how our results are affected if 
we make the extreme assumption of full kinetic equilib- 
rium in the bottom sector, i.e., set TZ = 1. We will 
again do this within the quasifree dissociation scenario 
including in- medium binding energies where the rates, 
and therefore the sensitivity, are largest. The results are 
summarized in Figs. 1201 and 1211 Even though the rela- 
tive increase in A^ q is substantial, its impact at RHIC 
(Fig. I20fl is still rather limited (~20%) since it becomes 
comparable to initial hard production (after nuclear ab- 
sorption) only after t ~ 2 fm/c when the temperature has 
dropped below 200 MeV and consequently the reaction 



rates have become small, cf. Fig- E] At LHC (Fig. 121(1 . 
even though the relative increase in is smaller than 
at RHIC (since the 6-quark relaxation time is smaller at 
larger temperatures), its absolute number is now closer to 
the actual T number, Ny(t), in the system, especially at 
intermediate times in the QGP evolution (t ~ 2-4 fm/c). 
Since the reaction rates at this stage are still sizable, 
the final number in the thermal-equilibrium scenario in- 
creases by about a factor of 2 over the off-equilibrium 
case. Nevertheless, the main feature of a strong suppres- 
sion is preserved (about a factor of 3), while the final 
yield is composed to about ^90% of regenerated T's. 

This should have noticeable consequences for p t -spectra 

1 1 i 1 1 1 I 1 1 1 i 1 1 1 i 1 1 1 I 1 1 1 i q 
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FIG. 21: (Color online) Thermal off-equilibrium effects at 
LHC for the suppression of T(ls) in central (6=1 fm) Pb- 
Pb collisions, using the quasifree process. The light set of 
curves indicates our results assuming thermalization of the 
open-bottom states, showing larger equilibrium abundances 
(dash-dotted curve) and a larger final T yield (solid line) than 
with thermal off-equilibrium effects (dark set of curves, cor- 
responding to Fig. Illl in the main text). 



(regenerated T's are expected to have much softer spec- 
tra than primordial ones, see, e.g., the recent analyses 
m Refs. ha, 0]), and presumabl y ev en more so for the 
elliptic flow (compare, e.g., Refs. [El] and indicating 
a difference of up to factor ~5 in V2 for the charmonium 
case). 
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